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the change in distribution function of the stellar mass production rate in individ- 



0> ■ ual galaxies with redshift. The evolution of the star formation rate distribution 

1^ • with redshift is interpreted in terms of the history of spiral galaxy formation, with 

I ' the disk component modelled as a single evolving entity, and the characteristic 

;_( ■ timescales, luminosities, and epochs varying according to galaxy type. The more 

^ ■ massive forming galaxies seen at 2; = 1 ^ 3 are identified as earlier type spirals, 

whose star formation rates are initially high and then decline rapidly at z < 1, 

i^j while for later type spirals and smaller mass irregulars the mass formation rates 

'j_j ■ at 2; < 1 are lower, and the formation process persists to redshifts much closer 



to the present epoch. We find that these models can be consistent with the data 
and fit well into a broad picture of other recent results if go=0.02 and many of 
the disks begin their growth at 2; <^ 3, but that they predict too many bright 
star formers at high z in fiat universes. 
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1. Introduction 



Disks of spiral galaxies are extremely fragile and must have formed primarily by con- 
version of gas to stars in an essentially invariant galaxy potential with only a small accretion 
of external stars (e.g., Toth & Ostriker 1992). Thus, the stellar history of this component 
may be approximately modelled as the evolution of a single entity, which greatly simplifies 
predicting the luminosity evolution. In the present Letter we trace the evolution of the 
distribution of star formation rates in galaxies with redshift to determine whether it is con- 
sistent with this type of pure luminosity evolution. The problem is relatively straightforward 
because beyond z ~ 0.3 the observed B band corresponds to ultraviolet light (A ^ 3500 A), 
which in all but the oldest galaxies is a measure of the ongoing rate of massive star for- 
mation. The B luminosity is then directly proportional to the rate of metal production in 
the individual galaxy, while the contribution of the ensemble of galaxies to the extragalactic 
background light measures the overall metal density production in the universe (Cowie 1988; 
Songaila, Cowie, & Lilly 1990; Lilly et al. 1996; Madau et al. 1996). Viewed this way, deep 
5-band redshift samples measure the distribution function of the stellar mass production 
rate in galaxies, which we will refer to as the M function, in contrast to infrared samples, 
which for all but the youngest galaxies measure the distribution of the mass in galaxies. 

In order to probe star formation rates (SFR) in individual galaxies to high redshifts, 
large and extremely deep 5-band redshift samples are required. Here we use a large, 95% 
complete sample of galaxies to -B = 24.5 with redshifts \,o z = 2.2, to measure the evolution 
of the M function with redshift. We find that smooth single-entity evolution models can 
work if go=0.02 and many disks begin to form only at later epochs, but predict too many 
bright star formers at high z in fiat universes. 



2. Data 

Cowie et al. (1996) described LRIS observations of a i? = 24.5 galaxy sample in the 
Hawaii fields SSA13 and SSA22 (26.2 arcmin^). However, even the more intensively studied 
SSA22 field was only 84% complete at 5 = 24-24.5, with 13 5 < 24.5 galaxies unidentified. 
Since extremely high completeness is required to determine the M function and the uniden- 
tified objects may correspond to higher redshift objects in the sample, we made intensive 
efforts to complete identifications of the remaining B < 24.5 galaxies in SSA22 using deeper 
optical (3500-10000 A) observations and near infrared spectroscopy. The improved S/N and 
additional spectral features yield robust identifications for all but four sample objects. 
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The SSA22 B < 24.5 sample is used with a nearly complete B < 24 sample in SSA13. 
In the B < 24 sample covering 26 arcmin^ there are 156 objects (124 galaxies) of which all 
but 4 are identified (97% complete), while in the B = 24 — 24.5 range covering 13 arcmin^ 
there are 47 objects (42 galaxies) of which all but 2 are identified (95% complete). The 
six unidentified galaxies have extremely blue colors but no strong [O ll] , with most likely 
redshifts 1.6 < z <2, where the [Oil] line has moved from the observed optical while many 
stronger UV features have not yet entered the blue wavelength range. However, some objects 
may lie at somewhat higher redshifts. 

The redshift-magnitude relation and median redshifts are shown in Figure |I|, along with 
the Autofib sample of Ellis et al. (1996). As at brighter magnitudes, the median B red- 
shifts, computed from the combination of the Autofib data plus the present sample, are 
well described by a model (solid line) in which the luminosity function remains invariant 
(Broadhurst, Ellis, & Shanks 1988; Colless et al. 1990; Cowie, Songaila, & Hu 1991; Glaze- 
brook et al. 1995a; Cowie et al. 1996), though this refiects a much more complex underlying 
situation. Dashed and dotted lines show expected magnitudes for a fiat F^ galaxy with 
Mab = -19.8 + 5 log/i65 for go=0.5 (dashed) and Mab = -20.5 + 5 log/igs for go=0.02 (dot- 
ted), which roughly map the upper envelope of the magnitude-redshift relation at 1 < z < 2. 
As we shall see below these correspond to star formation rates of roughly 9 h^^ M© yr^^ 
(go=0.5) and IGhg^ M© yr~^ (go=0.02). As Cowie et al. (1996) emphasized, there are few 
galaxies with such high star formation rates at 2; -C 1, and the envelope of the observed 
distribution lies well to the right of the curves at these redshifts, except for a few luminous 
AGN present in the samples. The brighter z = 2 — 3 galaxies lie to the left of the line, with 
mass formation rates up to a factor of 2.5 higher. (The Autofib and HDF data are not used 
in the remainder of the paper in order to avoid completeness problems.) 



3. Interpretation 

The present B < 24.5 sample provides two independent measurements of the massive 
star formation rate. The first is the rest-frame far UV (A ^ 3500 A) light, which for galaxies 
with significant ongoing star formation is a direct measure of the massive star formation 
rate, and hence of the metal production rate (Cowie 1988; Songaila et al. 1990). Translating 
this to a total star formation rate requires an uncertain assumption about the slope of the 
IMF. For a Salpeter IMF extending over the range 0.1 ^ 125 M© we may calibrate from the 
Bruzual & Chariot (1996) models to obtain, in the absence of internal extinction, 

M^100M,y.-( ^^^^^r°^).^_. ). (1) 

\ 6.7 X lO^y ergs s ^ Hz / 
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where a rest-frame wavelength of 2500 A is chosen to minimize extrapolations from the 
B band in lower redshift galaxies while also minimizing the slight contamination by older 
stars which is present at slightly longer wavelengths. Equation (|I]) may be derived from 
extremely simple physical arguments and is quite model independent except for the assumed 
IMF (Cowie 1988). 2500 A corresponds to B at z = 0.8 and I at z = 2.4, and may be 
interpolated at intermediate redshifts using approximate galaxy types inferred from the {B, 
I, K) galaxy colors, so that the 2500 A luminosity can be relatively accurately determined 
at least at the higher redshifts. The massive star formation rate can also be obtained 
from the line luminosities in the galaxies, which measure the production rate of ionizing 
photons. The most direct diagnostic is L(Ha), which can be measured in these galaxies 
from 2; = — > 2.3. For most of the galaxies, however. Ha is not measured, and we must 
draw on the secondary indicator [O ll] 3727 A, whose equivalent width is closely related to 
that of Ha by the relationship WA(Ha+[Nll]) = 2.3 Wa([0 ll]) in a wide variety of local and 
moderate redshift galaxies (Kennicutt 1992; Songaila et al. 1994). For local blue irregular 
galaxies (Gallagher, Bushouse, & Hunter 1989) 

M^100M,yr.(-M^). (2) 

yiO*'^ ergs s ^ ) 

but we derive an essentially identical calibration for the present galaxies by comparison with 
the more direct UV luminosity calibration, suggesting that this calibration is reasonable for 
most rapidly star-forming galaxies irrespective of mass and redshift. The [Oil] luminosity 
provides an independent check on the UV luminosity which is particularly valuable at the 
low redshift end (Cowie, Hu, & Songaila 1995). 

Internal extinction in the galaxies is a serious concern for both diagnostics, since escaping 
UV photons can be reduced significantly by this effect, while the internal production of Ha 
and [On] can be reduced if ionizing photons are extinguished. However, at least for the 
higher redshift blue galaxies there appears to be remarkably little extinction and this effect 
appears weak. The great bulk of the z > 0.8 galaxies are extremely fiat in the rest-frame 
2500 A— 21000 A color, with a median value of 0.9 in the AB system, and have high rest-frame 
[On] equivalent widths (with a median of 40 A). Since the intrinsic spectrum is expected to 
be fiat or falling with increasing frequency, this places very severe constraints on the amount 
of dust. For a Seaton (1979) law A2500 ~ ^(2500 A - 10000 A) ^ 7 E{B - V). Even if F^ 
is intrinsically fiat, A2500 = 0.9 on average, and at most half the light is lost to extinction. 
Thus, the maximum uncertainty introduced by not including an extinction correction is 
probably a factor of 2. To allow for the extinction we therefore assume that the luminosities 
in equations (P and (0) are reduced by a factor of 1.5 for a given mass formation rate. 



-5- 

The computed rest-frame 2500 A and [O ll] luminosities are shown versus redshift in Fig- 
ure 0. For go=0-02 and Hq = 65 km sec~^ Mpc~^ (Fig. 0a) the upper envelope of the 2500 A 
luminosity at z > 1 lies at approximately 10^^ h^^ ergs s~^ Hz~^ and just above 10^^ h^^ 
ergs s~^ for [Oil], which correspond to mass formation rates in the 15 — 30 h^^ Mq yr~^ 
range. No objects with much higher mass formation rates are seen. At lower redshifts the 
upper envelope falls rapidly. For go=0.5 we have used Hq = 50 km sec~^ Mpc~^ (Fig. ^b) to 
provide an acceptable age for the universe. With this smaller Hubble constant the maximum 
luminosities and mass formation rates are similar to those in the go=0.02 case. 

The distribution of M in various redshift bins, computed using the 1/V method, is shown 
for open and flat geometries in Figure ^ The star formation rate may be directly compared 
with the local star formation rate determined by Gallego et al. (1995) from their Ha survey, 
which covered redshifts z^O.045, for a 470 deg^ sample with EW(Ha+[Nll])> lOA. This is 
shown for z = 0.1 — 0.4 in the bottom right panel of Figure p. The data points lie roughly a 
factor of two higher than the Gallego et al. points, although both data sets are fully consistent 
in slope within the respective 2a errors and uncertainties. Interestingly, this same difference 
betwen the number density of local vs. distant galaxies is seen in the ii'-band count analysis 
of Huang et al. (1997), who suggested that there is a local deficiency of galaxies by a factor 
of 2 on scale sizes out to ~ 300 h~^ Mpc. Applying such a factor of 2 correction to the local 
SFR normalization would then directly match the data at z = 0.1 — 0.4. 

There is a significant evolution in the maximum mass formation rates with redshift for 
both open and fiat geometries, in the sense that galaxies with much higher formation rates 
are seen at the higher redshifts (e.g., Cowie et al. 1996). For z = 0.8 — 1.6 there are 27 
galaxies with M > 10 M0 yr~^ while none are seen a.t z < 0.8. Even allowing for the relative 
volumes the probability of their being drawn from the same distribution function is only 
6 X 10^^ for go=0.02 and 3 x 10^^ for go=0.5. However, the volume density of the rapid star 
formers is not high. For go=0.02 the volume density of objects with M > 10 h^^ Mq yr~^ 
is 5.4 X 10-^ hei Mpc^^ at 2 = 0.8 - 1.6 and 1.0 x 10"^ h^i Mpc'^ for z = l.Q - 3.2, while 
for go=0.5 and M > 10 h^o M© yr^-^ the densities are 6.2 x 10"^ h^i Mpc'^ at z = 0.8 - 1.6 
and 2.2 x 10"^ h^^ Mpc'^ for z = 1.6 ~ 3.2. 



4. Discussion 

In order to connect these results to the present epoch we need the mass distribution and 
universal mass density of the local galaxy sample, and these are best derived from the K- 
band observations. There is still some disagreement in the derived values of the local iiT-band 
luminosity function, and we adopt here the average of the samples of Mobasher, Sharpies, 
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& Ellis (1993), Glazebrook et al. (1995b), and Cowie et al. (1996), which corresponds to a 
Schechter function Mk* = -24.2 (+5 log^o /^es), « = -1-1, and 0* = 2x 10"^ hl^ Mpc"l The 
disagreement in Mk* is ~ 0.5 mags. Most of the light density of the universe is contained in 
objects near L*, with 50% of the light density in objects brighter than (L) = 0.7L*. Most of 
these brighter galaxies are earlier than Sb based on their B, I, K colors (Huang et al. 1997), 
and for consistency we have converted from K magnitudes to masses using the Bruzual & 
Chariot (1996) models, with F^, = 2.5 x 10^^ ergs s~^ Hz~^ at 21000 A corresponding to a 
stellar mass of 10^^ M0, whence L^: corresponds to 1.4 x 10^^ h^^ M©. The local stellar 
mass density of the universe is then 2.8 x 10^ h^^ Mq Mpc^'^. 

The present-day spectra of the spiral galaxies are well represented by evolutionary mod- 
els with exponentially decaying star formation rates. Bruzual & Chariot (1993) fit an Sb 
galaxy with a tsfr = 2 Gyr and an age of 8 Gyr, and an Sc galaxy with a tsfr = 7 Gyr 
and an age of 12 Gyr. In Figure ^ we show the predicted luminosities for such exponentially 
evolving galaxies, with masses corresponding to L* for tsfr = 1 Gyr and tsfr = 2 Gyr, 
and to 0.5 L^: for tsfr = 7 Gyr, for comparison with the observed luminosities. Because 
both L^ and the observed galaxy luminosities scale as h~^, the curves can only be adjusted 
via the x-axis conversion from redshift to time and by the choice of the geometry, so the 
absolute agreement in the normalization of observed and predicted luminosities with redshift 
shows impressive consistency. For go=0.02 we can fit the envelope with an Hq = 65 km see" 
Mpc~^ and a late formation epoch Zf = 3 model, but for go=0-5 we need a low Hubble 
constant (Ho~50 km sec~^ Mpc^^) and an early formation epoch (zj ~ 10). Basically, to 
match the fall of the maximum luminosity with redshift at z < 1, we require galaxy ages of 
very roughly 4 Gyr at z = 1. The expected luminosity envelope is better matched for the 
go = 0.02 case, while for go=0.5 there are fewer high luminosity galaxies at high redshift to 
form the more massive spiral galaxies. Later type galaxies, such as the 0.5 L* Sc shown by 
the dotted line in Figure 0, are also consistent with the observations provided they are not 
too massive, so that their luminosities are not too high at the lower redshift s. A requirement 
of the models is therefore that the later type galaxies should have lower K luminosities than 
the earlier types, with Scs being at least a magnitude fainter than Sbs, and Sds being con- 
siderably fainter still. The current if -band data are not yet adequate to robustly investigate 
type dependence in the luminosity function, but this prediction should shortly be testable. 

We can now quantify the problem with the go=0.5 geometry: Predicted luminosities 
of the early stages of near-L^, spirals are well above the detection limit (Fig. |^), at least 
for galaxies earlier than Sb, and these galaxies must have formed early {zf ^ 3) to provide 
large enough ages for their ultraviolet luminosities to match the observed redshift evolution 
and falloff at ^ ~ 1. Thus, in this case we should see a number density of galaxies in the 
z = 1.75 -^ 3 range which is comparable to the present-day number density (2.2 x 10~^ h^^ 
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Mpc"^ for L > 0.25 L* and 1.1 x 10"^ /igg Mpc"^ for L > 0.5 L*). Even if all the unidentified 
objects are allocated to this redshift range there are only 10 objects, corresponding to a 
number density of 4.3 x 10"'^ ± 1.5 x 10~^ /igg Mpc^^. For go=0.02 the number densities 
are similar, but the problem can be avoided because the formation can begin in this redshift 
range, so that number conservation is not required. In this case only a fraction of the disks 
can have begun to form at z > 2 with the remainder starting up at lower redshifts. 

In any given redshift interval we can compute the total universal mass density formation 
rate, p, in observed objects, which corresponds to a lower limit on p. Correcting for the 
incompleteness is not straightforward since it depends on the present-day type mix as a 
function of mass and the age at a given redshift. However, including only the observed 
objects and without correcting for faint end incompleteness, we find mass formation rates 
of M = (1.4 X 10-2, 1.7 X 10-^ 1.0 X 10-2,2.8 x lO^^) h^i Mq Mpc'^ yr-^ for z = (0.2, 
0.6, 1.2, 2.4) and go=0.02, though these numbers are not directly comparable because of the 
varying cutoff in M so that formation rates at high redshift are relatively underestimated. 
For objects with M > 10 M0 yr-^, where all objects would be detected at all of these 
redshifts, the average formation rates are (0.0, 0.0, 7.1 x 10-^,2.8 x 10-^) h^^ Mq Mpc-^ 
yr-^. The overall rates are sufficient to form the present-day mass, so that the open models 
do appear to provide a fully self-consistent description with the formation peaking at 2; ~ 1 
at least for the objects with higher star formation rates. 



5. Conclusion 

We conclude that in a go=0.02 universe we can understand much of the recent history 
of galaxy formation in terms of the pure luminosity evolution of various spiral galaxy types, 
provided only that Sc galaxies are at least a magnitude fainter in K than the average of 
galaxies earlier than Sb, with Sds yet fainter. Models with qQ=0.5 predict too many bright 
objects at higher redshifts. The conclusion that this works well for go=0.02 geometries but 
not for flat geometries is intimately connected with pure luminosity evolution modelling 
of the number counts and redshift distributions (e.g., Metcalfe et al. 1996), where go=0.5 
models tend to produce elbows in the counts and fail to reproduce the faint-end counts. The 
reason for this is that in the flat models the early stages of the spirals have near-constant 
and relatively bright magnitudes at all high redshifts (Cowie 1988) — the same problem we 
encounter here with a radically different approach. The go=0.02 models do not suffer from 
this problem and have much more latitude in available time and volume to resolve the high 
number count problem (Lilly, Cowie, & Gardner 1991), and also do achieve broad consistency 



(Metcalfe et al. 1996) with optical and IR number counts and redshift distributions in this 
geometry. 
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Fig. 1. — The redshift--B magnitude distribution for SSA22 B < 24.5 (solid boxes) and 
SSA13 B < 24 (diamonds). Unidentified objects are shown as crosses at nominal redshifts 
oi z = 2.6 (SSA22) and z = 2.7 (SSA13). The very small symbols show the large Autofib 
sample of Ellis et al. (1996). The few points lying outside of the well-populated regions are 
AGN. The large crosses show the median redshifts computed in half-magnitude intervals with 
±lo" error bars computed using the median sign method. (The Autofib data are included 
at i? < 23 only.) These are very well described by the no-luminosity-evolution model of 
Glazebrook et al. (1995a) (sohd line). The dashed (go=0.5) and dotted (go=0.02) lines show 
the i?-magnitude vs. redshift relation for a flat F^, galaxy which roughly matches the upper 
envelope of the observed objects in the z = 1 — 2 range. We also show (open triangles) the 
six B < 25, z > 1.6 galaxies which have been identified in the Hubble Deep Field to date 
(Cohen et al. 1996; Steidel et al. 1996b; Lowenthal et al. 1997). This rather heterogeneously 
selected sample shows rough consistency with the upper envelope of the B — z relation seen 
in the magnitude-selected sample, as do other galaxies selected with the ultraviolet break 
color techniques (Steidel et al. 1996a), though they may also contain a small number of more 
luminous galaxies. 
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Fig. 2. — The rest frame 2500 A luminosity (left panels) and [O II] luminosity (right panels) 
computed for two geometries — go=0.02 and Hq = 65 km sec^^ Mpc^^ (a) and go=0.5 and 
Ho = 50 km see" Mpc~^ (b). In the left-hand panels the dotted line shows the selection 
limit for a, B = 24.5 galaxy with a flat Fy spectrum. Unidentifled objects are shown as crosses 
at a nominal redshift of 2.5. Models are shown for exponentially decaying star formation 

Mq galaxy with r = 1 Gyr (solid line), and r = 2 



rates which would form a 1.4 x 10^^ h^^ 



Gyr (dashed line), or a 7 x 10"*^° h^^ Mq galaxy with r = 7 Gyr (dash-dot line). For the 
open geometry the galaxy evolution is started at 2;/ = 3 and for the flat geometry at Zf = 9. 
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Fig. 3. — The distribution function of the star formation rates as a function of redshift 
interval for go=0.02 (left panels) and go=0.5 (right panels). The vertical axis shows the 
number of galaxies per Mpc^ per logarythmic bins of 0.3. The errors are ±la based on the 
number of objects in each bin. The filled squares show the local {z^O.045) SFR of Gallego 
et al. (1995), plotted for the Salpeter IMF assumed for the present calculations. (The Ha 
luminosity corresponding to a star formation rate of 1 solar mass per year is a factor of 3 times 
higher than for the Scalo IMF, less rich in massive stars, used in the Gallego et al. paper.) 
Error bars (±lcr) for these points, based on the number of objects per bin, are comparable 
to the symbol size; for the endmost points, they are roughly twice the symbol size. The 
agreement in shape between these two curves is good, with the SFR for the z = 0.1 — 0.4 
range approximately a factor of two higher than for the Gallego et al. sample. Differences 
between the two curves are likely to arise due to differences in each sample's completeness 
in identifying star-forming galaxies, as well as a possible deficiency in local galaxies (Huang 
et al. 1997). 



